AFRL-VS-HA-TR-2006-1 004 

J.  Phys.  Chem.  A  2004,  108,  10965-10972  10965 


O  +  C„H2 n  Products  Detected  via  IR  Emission.  1.  O  +  C2H4 


D 

fcaaj 

o 

o 

o 


James  A.  Dodd,*  Eunsook  S.  Hwang, +  Karen  J.  Castle,*  and  Gary  D.  DeBoer5 

Air  Force  Research  Laboratory,  Space  Vehicles  Directorate, 

Hanscom  Air  Force  Base,  Massachusetts  01731-3010 

Received:  June  21,  2004;  In  Final  Form:  August  31,  2004 

Collisions  between  oxygen  atoms  and  ethene  have  been  investigated  by  using  infrared  emission  detection  of 
the  chemiluminescent  product  species.  A  prototypical  alkene,  ethene,  nonetheless  exhibits  numerous  reaction 
pathways  in  reactions  with  O  atoms.  Oxygen  atoms  were  formed  via  photolysis  of  SO2  in  the  presence  of 
C2H4,  and  the  resultant  IR  emissions  in  die  900— 3000-cm"1  spectral  region  were  detected  by  using  a  time- 
resolved,  step-scan  Fourier  transform  spectrometer.  A  Welsh  cell  mirror  arrangement  was  used  to  maximize 
the  collection  efficiency  of  the  product  IR  emissions.  Vibrationally  excited  products  such  as  CO,  CO2,  HCO, 
and  H2CO  have  been  identified,  with  CO  and  C02  being  the  dominant  IR  emitters.  The  time-evolving  CO 
and  CO2  spectra  have  been  characterized  with  respect  to  the  SO2  and  C2H4  partial  pressures  and  laser  fluence. 
The  rate  constants  for  vibrational  relaxation  of  CO2  high-t/  population  by  C„H2„  (n  =  2-4)  are  in  the  mid- 
10“ 12  cm3  s-1  range;  SO2  is  a  very  inefficient  relaxer.  A  chemical  kinetics  code  has  been  used  to  model  the 
chemistry  and  identify  the  operative  reaction  mechanisms,  including  the  effects  of  secondary  chemistry. 


1.  Introduction 

Hydrocarbons  are  not  present  to  any  significant  extent  in  the 
Earth’s  upper  atmosphere;  those  produced  through  anthropo¬ 
genic  sources  on  the  surface  are  broken  down  by  solar  radiation 
at  much  lower  altitudes.  However,  hydrocarbons  can  be 
introduced  into  the  upper  atmosphere  in  situ  in  the  vicinity  of 
a  space  vehicle,  for  instance  through  surface  outgassing  or 
leakage  of  unburned  fuel.  For  a  vehicle  in  low-earth  orbit  with 
a  velocity  of  7.8  km  s'1,  the  resultant  collisions  with  ambient 
atmospheric  species  take  place  with  center-of-mass  collision 
energies  Ecom  on  the  order  of  several  hundred  kJ  mol-1.  For  a 
suborbital  vehicle  with  an  apogee  in  the  100—150  km  range, 
the  interaction  energies  are  much  lower,  in  the  10—50  kJ  mol-1 
range,  though  still  significantly  superthermal. 

Above  about  80-km  altitude,  ground-state  0(3F)  atoms 
become  increasingly  prevalent  owing  to  solar  VUV  photodis¬ 
sociation  of  O2.  The  O-atom  density  peaks  at  about  100  km, 
with  concentrations  exceeding  10u  cm-3,1  and  exceeds  the  O2 
density  for  altitudes  above  about  120  km.  The  presence  of  O 
atoms  is  important  in  part  because  O  is  much  more  reactive 
than  the  other  prevalent  neutral  constituents,  N2  and  O2.  The 
high  collision  energies  that  are  possible— for  instance,  310  kJ 
mol-1  for  O— C2H4  collisions  at  the  orbital  velocity— serve  to 
stimulate  two  phenomena:  first,  nonreactive  excitation  of  the 
internal  modes  of  the  hydrocarbon,  and  second,  chemical 
reaction  by  surmounting  the  potential  barrier  governing  the 
reaction  coordinate,  resulting  in  various  molecular  and  radical 
products.  Lower  energy  processes,  typically  involving  secondary 
chemistry,  can  also  occur  in  the  vehicle’s  wake,  downstream 
of  the  initial  interaction  region  near  the  ram  surface.  These 
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downstream  radical— radical  and  molecule— radical  reactions 
often  occur  at  gas  kinetic  rates  with  little  or  no  temperature 
dependence,  giving  rise  to  fast  processes  even  in  the  relatively 
cool  wake  region.  In  both  cases,  the  resultant  physical  excitation 
and  chemical  reaction  processes  yield  rovibrationally  excited 
products  that  often  radiate  in  the  infrared.  The  structure  of  this 
emission  could  provide  a  signature  as  to  the  type  and/or  origin 
of  the  vehicle  in  question. 

In  this  report,  we  provide  results  from  a  recently  initiated 
experiment  detecting  the  products  of  C2H4— O  collisions.  The 
C2H4-O  system  was  chosen  for  several  reasons:  (1)  alkenes 
are  common  components  of  missile  fuels,  (2)  the  chemistry, 
even  for  the  simplest  alkene,  is  complex  and  results  in  a  host 
of  IR-active  products,  and  (3)  the  study  complements  work  being 
performed  in  our  laboratory  with  laser-induced  fluorescence 
detection  of  the  product  HCO  radical.2  Two  long-term  goals  of 
this  work  are  to  determine  detailed  product  state  information 
and  to  observe  the  role  of  translational  energy  in  the  reactions. 
Reactions  of  alkenes  with  O  atoms  have  undergone  a  large 
amount  of  study  within  both  the  combustion3  and  chemical 
dynamics  communities.4,5  Figure  1  shows  a  schematic  of  the 
principal  reaction  pathways  available  to  O  +  C2H4  encounters. 
The  main  product  channels  are  an  H-atom— vinoxy  radical  pair 
(H  +  CH2CHO)  and  a  formyl— methyl  radical  pair  (HCO  + 
CH3).  There  is  some  disagreement  as  to  the  partitioning  between 
these  channels,4,6-8  though  each  accounts  for  roughly  half  of 
the  overall  product  yield.  According  to  at  least  one  study,  the 
methylene— formaldehyde  (CH2  +  CH2O)  channel  accounts  for 
a  few  percent  of  the  yield.7  The  further  unimolecular  dissociation 
of  CH2CHO  to  several  different  products  has  also  been 
studied,6,9  as  has  that  of  HCO  to  H  +  CO,  both  experimen¬ 
tally10-12  and  theoretically.13-16  The  H— CO  dissociation  channel 
is  shown  with  its  approximately  96  kJ  mol-1  barrier,  together 
with  the  lowest  four  product  CO  vibrational  levels. 

In  this  work,  vibrationally  excited  products  such  as  CO,  CO2, 
HCO,  and  H2CO  have  been  detected  via  emission,  using  the 
technique  of  time-resolved  Fourier  transform  infrared  (FTIR) 
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Figure  1.  Reaction  path  schematic  of  the  O  +  C2H4  reaction,  including 
the  6.65  kJ  mol-1  room  temperature  activation  barrier  (ref  27).  The 
principal  direct  reaction  channels  are  H  +  CHiCHO  and  HCO  +  CH3, 
where  each  represents  about  half  of  the  overall  product  yield.  The  H2 
+  CH20  channel  accounts  for  a  few  percent  of  the  yield  (ref  7).  The 
H-CO  dissociation  channel  is  shown  with  its  approximately  96  kJ 
mol-1  barrier,  as  well  as  the  lowest  four  CO  vibrational  level  energies. 

spectroscopy.  While  FTIR  emission  detection  is  a  versatile  and 
powerful  means  of  observing  the  IR-active  products  in  a  reaction 
mixture,  it  has  the  disadvantage  of  low  sensitivity  relative  to 
other  detection  methods.  Thus,  the  species  partial  pressures  are 
necessarily  relatively  high  to  ensure  adequate  signal-to-noise 
ratios  for  detection  purposes.  Under  these  conditions,  secondary 
chemistry  plays  a  significant  role  in  the  production  of  IR-active 
species.  To  better  understand  the  chemistry  occurring  in  the 
reaction  cell  following  the  excitation  pulse,  a  chemical  kinetics 
model  was  employed  to  predict  species  concentrations  as  a 
function  of  time  after  laser  excitation,  including  the  effects  of 
secondary  chemistry  following  the  initial  O  +  C2H4  reaction. 

2.  Experimental  Section 

Oxygen  atoms  were  generated  through  the  193-nm  photolysis 
of  SO2  with  use  of  95-mJ  light  pulses  from  a  30-Hz  Lambda- 
Physik  Lextra  200  excimer  laser.  The  S02  photolysis  produces 
a  relatively  peaked  distribution  of  O-atom  energies,  centered  at 
27.0  kJ  mol-1  laboratory  frame  energy,17-20  or  17.2  kJ  mol-1 
COM  energy  for  the  O-C2H4  interaction.  Mass  flows  of  15 
standard  atm  cm3  min-1  (seem)  SO2,  15  seem  C2H4,  and  100 
seem  Ar  bath  gas  were  passed  through  a  reaction  cell  at  a  total 
pressure  of  500  mTorr.  The  cell  was  fitted  with  UV-fused  silica 
windows  for  the  laser  excitation  pulses  and  a  CaF2  window, 
positioned  at  90°,  for  sampling  the  IR  chemiluminescence.  The 
gas  mix  was  introduced  and  evacuated  via  opposing  2.54-cm- 
diameter  pipes  separated  by  10  cm.  The  cell  contained  a  pair 
of  10.2-cm-diameter,  10.2-cm  radius-of-curvature  aluminum- 
coated  mirrors  arranged  in  the  Welch  cell  configuration21-24 
for  collection  of  IR  chemiluminescent  emission.  A  portion  of 
the  Ar  bath  gas  was  used  to  purge  the  sidearms  supporting  the 
laser  windows  and  the  mirrors. 

The  IR  emission  exiting  the  cell  reflected  off  two  aluminum 
mirrors  into  a  Bruker  DFS/66  step-scan  spectrometer,  capable 
of  obtaining  simultaneous  spectrally  and  temporally  resolved 
data.  For  these  experiments  the  nominal  spectral  resolution  was 
either  2  or  15  cm-1,  with  250  time  samples  of  5-/us  time 
resolution  collected  following  laser  excitation.  Emission  data 
from  up  to  100  laser  shots,  or  3.3-s  real  time  at  the  30-Hz 


TABLE  1:  IR  Emissions  in  the  1000—3500  cm-1  Spectral 
Range,  from  193-nm  Laser  Excitation  of  Sample  Gas 
Mixtures 


principal  IR 

gas(es)  present"  emissions,  cm-1 


n-GiHg 

S02 

SO2  +  CO 
S02  +  11-C4H8 


2750-3350 
1150,  1360 
1150, 1360 
1150, 1360 
1900-2300 
2400-3000 


emitting  bands 


aliphatic  C-H  stretch  modes 
SO(y),  S02(vi,v3) 

SO(i>),  S02(vi,v3) 

SO(y),  S02  (vi,v3) 

CO(y),  C02(v3) 

HCO(v,),  H2CO(v,,v5) 


"  In  addition  to  Ar  bath  gas. 


repetition  rate,  were  averaged  at  each  mirror  position.  A  1-mm 
square  HgCdTe  detector  was  employed  as  either  “openband” 
to  cover  the  900—3500  cm-1  range  or  with  a  cold  band-pass 
filter  limiting  transmission  to  the  1950—2550  cm-1  range 
(“filter”  data).  The  real  time  required  to  obtain  a  data  matrix  of 
two-dimensional  intensity  information  is  proportional  to  the 
wavenumber  band-pass  divided  by  the  spectral  resolution.  Thus, 
acquiring  a  filter  data  matrix  could  be  accomplished  in  about 
one-quarter  the  time  required  for  an  openband  data  matrix  with 
equivalent  S/N. 

With  the  3.3  L  s-1  pumping  speed,  the  e-fold  pumpdown 
time  in  the  5.6-L-volume  reaction  cell  is  equal  to  several 
seconds.  However,  the  volume  in  the  cylinder  between  the  input 
and  output  pipes,  50  mL,  is  evacuated  on  a  time  scale  of  tens 
of  milliseconds.  Since  the  reactant  gases  are  introduced  wholly 
through  the  input  pipe,  while  the  Ar  bath  gas  is  introduced 
mainly  through  the  sidearms  and  more  remote  portions  of  the 
cell,  the  latter  time  scale  is  likely  to  be  a  more  accurate  measure 
of  the  residence  time  of  the  reactive  gas  mixture.  A  buildup  of 
certain  reaction  products  could  complicate  the  chemistry,  for 
instance  through  the  photolysis  of  CH2O.  To  test  whether  this 
was  affecting  the  present  results,  the  experiment  was  run  at  10 
Hz  in  addition  to  the  normal  30  Hz;  the  slower  repetition  rate 
significantly  decreases  any  buildup  of  reaction  products  between 
laser  pulses.  Decreasing  the  repetition  rate  caused  no  discernible 
effect  in  the  observed  time-dependent  IR  spectra  for  the  gas 
flows  used  herein. 

Table  1  summarizes  a  series  of  control  experiments  to  check 
for  the  source  of  the  emissions  ascribed  to  chemical  reaction 
between  O  and  C2H4.  All  the  species  were  dilute  in  Ar  bath 
gas  in  a  roughly  1 : 1  ratio  when  more  than  one  minor  species 
was  present.  The  alkene  n-but-l-ene  (n-C<|Hg)  used  for  these 
experiments  has  a  much  larger  193-nm  absorption  coefficient 
than  C2H4,  with  the  absorption  cross  section  <7i93(n-C4Hg)  = 
5.9  x  10-19  cm2  25  equal  to  about  5%  of  amfSCh)  =  1.2  x 
10-17  cm2.  Thus,  the  photolysis  of  neat  n-C^g  provides  a 
stringent  test  regarding  the  possible  formation  of  chemilumi¬ 
nescent  hydrocarbon  fragmentation  and/or  reaction  products.  As 
seen  in  Table  1,  193-nm  photolysis  of  dilute  n-CaHg  produced 
only  short-lived  emission  (<40  /ts)  in  the  2750—3350  cm-1 
region,  corresponding  to  aliphatic  C-H  stretch  emission  at 
higher  wavenumber  values  than  those  relevant  to  the  analysis 
below.  Photolysis  of  neat  S02  or  SO2— CO  mixtures  gave  rise 
to  SO2  photoproducts  only.  Finally,  photolysis  of  S02— n-C+Hg 
mixtures  produced  the  CO  and  CO2  emissions  in  the  1900— 
2300  cm-1  region  analyzed  herein. 

3.  Results  and  Discussion 

3.1.  IR  Radiance  Observations.  Figure  2  provides  a  three- 
dimensional  “waterfall”  view  of  openband  emission  spectra  at 
selected  time  delays  following  laser  excitation  of  the  SO2— 
C2H4— Ar  gas  mixture.  Various  IR  bands  are  observed,  consis- 
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Figure  2.  Time  evolution  of  openband  IR  emission  spectra  resulting 
from  193-nm  laser  excitation  of  an  SC^-CoF^-Ar  gas  mixture,  in  the 
1000-3000  cm-1  wavenumber  region  with  15-cnT1  resolution.  Sample 
spectra  are  shown  for  the  5-640  /is  delay  time  range  on  a  log  scale, 
with  each  successive  spectrum  delayed  an  additional  factor  of  2  in  time. 
The  spectra  have  been  corrected  for  detector  response.  See  text  for 
description  of  features. 
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Figure  3.  Time  evolution  of  openband  IR  emission  spectra  resulting 
from  193-nm  laser  excitation  of  an  SCh-CzKt-Ar  gas  mixture,  in  the 
1800-2400  cm-1  wavenumber  region  with  15-cnr1  resolution.  Fun¬ 
damental,  hot  band,  and  combination  v3  emission  are  observed  from 
vibrationally  excited  CO2,  together  with  overlying  fundamental  and  hot 
band  CO  (y  —  v  —  1)  emission.  The  CO2  emission  is  partially  absorbed 
by  about  1  m  of  room  air,  especially  for  wavenumber  values  larger 
than  2300  cm-1.  The  data  have  been  corrected  for  detector  response. 

tent  with  emission  from  the  chemiluminescent  CO  (v  —*  v  —  1) 
and  CO2  v3  (1900—2300  cm-1),  H-CO  v\  (2438  cm-1),  and 
H2CO  Vi  and  V5  (2450—3000  cm-1)  modes,  as  well  as  from 
the  precursor  species  SO2  vi  (1151  cm-1)  and  v3  (1362  cm-1) 
and  SO  (v  — *  v  —  1)  (1150  cm-1)  modes,  the  last  of  which  is 
overlapped  with  SO2  v\.  There  is  also  evidence  of  parent  C2H4 
v1  (949  cm-1)  and  V12  (1444  cm-1)  emission.  A  disadvantage 
of  the  FTIR  detection  technique,  at  least  at  low  resolution,  is 
that  intense  emissions  can  obscure  others  that  are  not  as  intense. 
For  example,  the  very  strong  CO  and  CO2  features  could  be 
concealing  weak  HCO  v2  emission  at  1868  cm-1. 

Figure  3  shows  details  of  the  spectral  region  dominated  by 
C02  (y3  — *  vs  —  1)  asymmetric  stretch  and  CO  (v  -*  v  —  1) 
fundamental  emissions,  using  a  color-coded  contour  plot 
representation  of  the  openband  spectrum.  Emission  from  highly 
vibrationally  excited  CO2  and,  to  a  lesser  extent,  CO  results  in 


Wavenumber  (cm*1) 

Figure  4.  Top:  Emission  spectrum  resulting  from  193-nm  laser 
excitation  of  an  SO2— C2H4— Ar  gas  mixture,  in  the  1925-2350  cm-1 
wavenumber  region  with  2-cm-1  resolution.  The  spectrum  was  obtained 
100  after  the  laser  pulse.  Bottom:  Synthetic  spectral  fit  consisting 
of  sharp  CO  lines  superimposed  on  a  broad  CO2  emission  feature.  The 
bold  line  highlights  the  C02  contribution.  See  text  for  discussion. 

hot  and  combination  bands  that  shift  the  radiance  to  lower 
wavenumber  values.  The  CO  v  populations  are  collisionally 
relaxed  relatively  slowly,  as  evidenced  by  the  “fingers”  at  2100 
and  2175  cm-1  extending  to  long  delay  times.  These  are  the  P 
and  R  branches  of  the  v  =  1— *0  fundamental  band,  which  is 
the  last  band  to  be  relaxed  as  the  vibrational  populations  cascade 
downward.  Preliminary  synthetic  spectral  modeling  of  the  CO 
v  emissions  is  consistent  with  a  Boltzmann-weighted  vibrational 
distribution,  with  an  effective  vibrational  temperature  of  about 
4000  K  at  small  delay  times.26 

For  the  C02  emission  profile,  preliminary  modeling  is  also 
consistent  with  a  highly  excited  vibrational  distribution  closely 
following  the  excitation  laser  pulse.26  At  early  times  the  emitting 
populations  have  undergone  negligible  vibrational  relaxation, 
and  the  radiance  in  a  spectral  region  that  preferentially  samples 
CO2  can  be  taken  as  proportional  to  the  C02  production  rate. 
For  the  most  highly  vibrationally  excited  CO2  populations, 
corresponding  to  those  radiating  in  the  longest  wavelength 
portion  of  the  spectrum,  there  is  negligible  feed  from  higher 
lying  populations.  This  lack  of  feed  simplifies  the  vibrational 
relaxation  measurements  discussed  below.  The  1900—2000 
cm-1  spectral  region  is  convenient  in  monitoring  the  CO2 
population  evolution:  the  emission  is  strong  enough  to  easily 
detect,  it  is  occurring  from  high-lying  vibrational  levels  of  C02, 
and  radiance  from  CO2  dominates,  especially  at  early  times. 
Spectral  fitting  indicates  that  in  the  1900—2000  cm-1  region, 
the  CO  v  radiance  contributes  at  most  a  few  percent  over  the 
first  400  /is.  At  later  times  the  contribution  by  radiating  CO  v 
increases,  but  the  overall  intensity  becomes  very  small  and  the 
decay  has  negligible  effect  on  the  results  of  the  CO2  population 
kinetic  analysis.  Figure  4  shows  a  sample  synthetic  spectral  fit 
to  a  100-/rs  delay  spectrum  for  an  excited  SO2— C2H4— Ar  gas 
mixture.  The  model  spectrum  consists  of  predicted  emission 
from  CO  vibrational  populations  characterized  by  a  2500  K 
effective  temperature,  and  from  a  distribution  of  CO2  rovibra- 
tional  populations  with  energies  up  to  40  000  cm-1.  As  is 
evident,  virtually  all  of  the  radiance  for  wavenumber  values 
<2000  cm-1  is  due  to  emission  from  excited  C02. 

3.2.  O-Atom  Thermalization  and  Chemistry.  The  COM 
energy  of  the  O  +  C2H4  system  afforded  by  photolysis  of  SO2 
is  equal  to  17.2  kJ  mol-1,  which  is  sufficient  to  overcome  the 
6.65  kJ  mol-1  room  temperature  activation  barrier  (Figure  l).27 
The  O— C2H4  reaction  rate  constant  is  strongly  temperature 
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TABLE  2:  Rate  Constants  Employed  in  the  Chemical  Kinetics  Model,  in  Units  of  cm3  s  1 


C2H4  +  O  —  CH2CHO  +  H 

3.75  x  10~'3 

CH3  +  CH3CO  —  CH3COCH3 

7.01  x  lO-" 

C2H4  +  O  — CH3  +  HCO 

3.30  x  10-'3 

CH3  +  CH3CO-  CH4  +  CH2CO 

1.01  x  lO"11 

c2h4  +  O  —  CH2CO  +  h2 

7.50  x  10"15 

CH3  +  CH2CO-*CO  +  C2H5 

8.30  x  10-'2 

C2H4  +  O  —  CH>  +  CH20 

4.50  x  10-14 

CH2CO  +  CH2  —  CH3  +  HCO 

1.00  x  10-'4 

C2H4  +  O  —  C2H3  +  OH 

2.56  x  10-,s 

CH3CO  +  H  —  CH3  +  HCO 

3.56  x  10-" 

CH2O  +  O  -  HCO  +  OH 

1.73  x  10-13 

CH3CO  +  H  —  CH2CO  +  H2 

1.92  x  10"" 

CH2O  +  OH-HCO  +  H2O 

1.00  x  10-" 

CH3CO  +  OH  —  CH2CO  +  H20 

2.01  x  10-" 

CH3CO  +  O  —  CH2CO  +  OH 

6.40  x  10-" 

CH3CO  +  o2  —  ch3o  +  co2 

4.51  x  10-’3 

CH3CO  +  O  —  CO  +  ch3 

2.56  x  lO-10 

CH2CO  +  H  —  CO  +  CH3 

7.32  x  lO-13 

CH3CO  +  CHi—  CH2CO  +  CHj 

3.01  x  10-" 

HCO  +  H  —  CO  4-  H2 

1.13  x  10-'° 

CH2  +  O-HCO  +  H 

5.01  x  10-" 

HCO  +  OH  —  CO  +  H20 

1.69  x  lO- 10 

CH2  4-  O  —  CO  +  H2 

8.00  x  10-" 

HCO  4-  CH2  —  CO  +  CH3 

3.01  x  lO-" 

CH2  +  O-CO  +  2H 

1.20  x  10-10 

HCO  4-  CH3CO  —  CHjCOCHO 

3.01  x  lO-" 

CH2  +  O  —  CH  +  OH 

7.20  x  10“12 

HCO  +  HCO  —  CH20  +  CO 

5.00  x  10~" 

CH3  +  O-CH2O  +  H 

9.50  x  10-" 

HCO  4-  02  —  CO  4-  H02 

3.86  x  lO'12 

CH3  +  O  —  CO  4-  H2  +  H 

4.00  x  lO-" 

HCO  +  HCO  —  2CO  +  H2 

3.64  x  lO-" 

CH3  +  O  —  CH30 

2.60  x  lO"14 

HCO  +  HCO  —  CHOCHO 

5.00  x  lO-" 

HCO  +  O  — •  CO  +  OH 

5.00  x  lO-11 

CH3O  +  CH3CO  —  CH20  4-  CH3CHO 

1.00  x  lO-" 

CH3CHO  +  O  —  CH3CO  +  OH 

3.90  x  lO"13 

CH3O  +  CH3CO  —  CH3OH  +  CH2CO 

1.00  x  lO-" 

HCO  +  O-CO2  +  H 

5.00  x  10-" 

CH2OCH2  —  CH3  +  HCO 

1.00  x  10->4 

OH  +  O  — *  O2  +  H 

2.50  x  10-" 

CH3CHO  4-  CH3CO  —  (CH3)2CO  +  HCO 

2.84  x  10-13 

OH  +  CO-CO2  +  H 

1.26  x  10-'3 

CH3O  4-  HCO  —  CH3OH  4-  CO 

1.50  x  10->° 

CH2CO  +  CH2  —  C2H4  +  CO 

2.09  x  10- 10 

H  +  HCCO  —  CO  +  CH2 

2.49  x  lO-10 

CH3  +  OH  — ■  H2  +  HOCH 

6.13  x  lO”13 

H  4-  CH2  —  H2  4-  CH 

4.92  x  lO"10 

CH3  +  OH  —  CH2OH  +  H 

1.31  x  10-" 

CH3  +  CH2OH  —  CH20  4-  CH4 

1.41  x  lO-10 

CH3  +  OH  —  CH3OH 

1.69  x  10-‘° 

CH3  +  S02  —  CH3SO2 

2.91  x  lO'13 

CH3  +  OH  -  H20  +  CH2 

1.13  x  lO"12 

SO  +  OH  —  SO2  +  H 

8.60  x  10-" 

CH3  +  CH2  —  C2H4  +  H 

7.01  x  lO-" 

CH2CO  +  OH  —  CH20  +  HCO 

4.65  x  10-" 

CH3  +  CH2OH  —  C2H5OH 

2.01  x  10-" 

02  +  SO  —  so2  +  O 

8.65  x  lO-17 

CH3  +  CHjCO  —  C2H6  +  CO 

5.40  x  lO-" 

CH2CO  +  OH  -*  H20  4-  HCCO 

1.19  x  lO'13 

CH3  +  HCO—CH3CHO 

3.01  x  10-" 

CH2CO  +  OH  —  CH2OH  +  CO 

7.10  x  lO-12 

CH3  +  HCO  —  CH4  +  CO 

4.40  x  10-" 

dependent,28  with  a  room  temperature  reaction  rate  of  7.5  x 
1CT13  cm3  s'1.  Prior  to  any  collisional  relaxation  of  the  nascent 
O  atoms,  the  O-C2H4  interaction  energy  corresponds  to  a 
kinetic  temperature  of  about  2000  K,  for  which  the  reaction 
rate  constant  is  45  times  higher,  or  3.4  x  10-n  cm3  s-1.  A 
small  portion  of  the  reaction  products  occurs  through  collisions 
involving  fast  O  atoms.  Ethene  is  present  in  a  0. 1 1  mol  fraction 
in  the  gas  mixture  and  thus  about  1%  of  the  initially  formed 
fast  O  reacts  prior  to  thermalization.  It  can  be  seen  from  Figure 
2  that  aside  from  prompt  emission  from  SO  v  and/or  SO2  Vi 
near  1150  cm-1,  the  majority  of  the  emission  has  a  delayed 
appearance,  consistent  with  primary  chemistry  induced  by 
thermalized  O  atoms  and/or  secondary  chemistry.  The  rate  of 
thermalization  of  the  energetic  O  atoms  can  be  estimated  by 
using  the  data  of  Matsumi  et  al.,29  who  found  that  41-kJ  mol-1 
O CD)  atoms  are  translationally  thermalized  by  Ar  with  a  near¬ 
gas  kinetic  rate  constant  of  1.28  x  10-10  cm3  s-1.  Assuming 
0(3P)  is  thermalized  by  Ar  with  a  similar  efficiency,  and  that 
Ar  is  representative  of  the  gas  mixture  in  general,  at  0.5-Torr 
pressure  the  e-fold  time  for  thermalization  is  about  0.5  /us.  Thus 
the  IR  emissions  from  the  chemical  reaction  products  discussed 
herein  arise  principally  from  thermal  chemistry. 

The  reactions  of  O  atoms  with  small-chain  alkenes  have  been 
studied  extensively  in  recent  years  by  Bersohn  and  co- 
workers.9,30-33  The  mechanism  involves  addition  of  0(3P)  to 
the  alkene  double  bond  to  form  a  triplet  biradical.  The  biradical 
can  either  dissociate  directly  to  produce  vinoxy  radical  and  an 
H  atom  or  intersystem  cross  to  the  singlet  state.  In  the  singlet 
state  the  barrier  to  H-atom  transfer  is  lowered  dramatically,  and 
the  energized  intermediate  can  give  rise,  through  a  combination 
of  H-atom  migration  and  dissociation,  to  a  host  of  daughter 
species. 

According  to  the  literature,  CO  and  C02  are  not  nascent 
products  of  the  O  +  C2H4  reaction.  However,  both  of  these 
species  can  be  formed  through  highly  exothermic  secondary 
reactions,  with  AHnn  values  in  the  —290  to  —470  kJ  mol-1 


range.  The  principal  formation  processes  as  identified  by  the 
kinetic  model  described  below  are  given  by 

O  +  CH3  —  CO  +  H2  +  H  (AHm  =  -287.4  kJ  mol"1) 

O  +  HCO  — CO  +  OH  (AHm  =  -363.8  kJ  mol"1)  (2) 
H  +  HCO  — -  CO  +  H2  (AHna  =  —371.6  kJ  mol-1)  (3) 
0  +  HC0-*C02  +  H  (A/7ran  =  -467.8  kJ  mol-1)  (4) 


These  energy  releases  are  sufficient  to  populate  highly  excited 
vibrational  levels  of  the  product  CO  and  C02  even  in  the  absence 
of  translationally  excited  reactants.  The  details  of  the  CO  and 
CO2  vibrational  population  distributions  and  formation  mech¬ 
anisms  will  be  discussed  in  a  forthcoming  publication.26 

3.3.  Modeling  of  Chemistry.  To  identify  the  chemical 
formation  mechanisms  of  the  principal  IR-emitting  products, 
including  CO  and  CO2,  a  chemical  kinetics  model  was 
developed  that  employs  the  commercial  REACT  computer 
code.34  Given  initial  concentrations  and  standard  chemical 
kinetic  inputs,  the  model  integrates  the  coupled  differential 
equations  describing  the  set  of  chemical  reactions  and  predicts 
the  species  concentrations  as  a  function  of  time,  with  an 
integration  step  size  of  5  /is.  All  reactions  were  characterized 
by  their  295  K  rate  constants.  A  total  of  40  species  and  65 
reactions  were  included  in  the  model  with  the  parameters 
obtained  from  multiple  sources,  notably  from  the  review  of 
Baulch  et  al.28  Table  2  details  the  chemical  reactions  and 
associated  rate  constants  used  in  the  model. 

Figure  5  shows  the  predicted  time-dependent  species  con¬ 
centrations  for  the  gas  mixture  described  in  the  Experimental 
Section,  assuming  SO2  photolysis  takes  place  at  t  —  0,  initiating 
the  chemistry.  About  a  third  of  the  SO2  is  dissociated  by  the 
typical  laser  fluence  used.  The  model  can  easily  predict  the 
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Figure  5.  Species  concentrations  predicted  by  the  chemical  kinetics 
model.  The  reaction  chemistry  is  initiated  at  time  t  =  0  through  the 
193-nm  photolysis  of  1.9  x  1015  cm-3  of  SO2  to  produce  SO  +  Of3/5), 
with  the  O  atoms  reacting  with  1.9  x  10ls  cm"3  C2H4  to  produce 
numerous  products.  Approximately  40  chemical  species  and  65  separate 
reactions  were  included  in  the  model.  Reaction  rates  and  product 
branching  ratios  were  taken  principally  from  ref  28.  The  model  assumes 
thermal  chemistry  throughout. 

product  species  formation  rates  as  a  function  of  the  inputted 
initial  O  and  ethene  concentrations.  These  predictions  are 
compared  with  the  experimental  observations  in  the  discussion 
below.  In  this  study  we  focus  on  CO2  production,  probed  by 
emission  from  highly  excited  vibrational  levels  in  the  1900— 
2000  cm-1  spectral  region. 

3.3.1.  Dependence  of  CO 2  Emission  Intensity  on  [OJinMai-  The 
band  intensity  for  excited  CO2,  integrated  between  1900  and 
2000  cm-1  and  0—50  ps,  is  shown  as  a  function  of  O-atom 
concentration  on  a  log— log  scale  in  Figure  6a  for  the  density 
range  appropriate  for  the  experiment.  Two  experimentally 
determined  curves  are  shown,  varying  the  laser  fluence  and  the 
SO2  concentration,  since  both  serve  to  increase  the  initial  O-atom 
concentration.  The  integrated  radiance  reaches  a  maximum  at 
a  delay  time  that  is  independent  of  fluence  or  SO2  concentration, 
and  thus  the  time  integration  is  a  good  representation  of  the 
rate  of  formation  of  vibrationally  excited  CO2  product.  The  slope 
for  the  log— log  plot  provides  an  estimate  for  the  order  of  the 
process  that  gives  rise  to  the  IR-emitting  CO2.  In  Figure  6b, 
the  model  prediction  for  CO2  chemical  production  as  a  function 
of  the  initial  O-atom  density  is  shown.  The  average  slope 
predicted  by  the  chemical  model,  1.42,  falls  midway  between 
the  experimental  values  of  1.31  and  1.61,  confirming  that  the 
dependence  is  intermediate  between  linear  and  quadratic.  The 
cause  of  the  discrepancy  between  the  experimental  laser  fluence 
and  SO2  concentration  results  is  not  clear,  though  the  slopes 
marginally  agree  to  within  the  error  bar  limits. 

The  formation  of  CO2  is  dominated  by  the  following  two- 
step  mechanism,  shown  with  the  room  temperature  rate  con¬ 
stants: 

O  +  C2H4  —  HCO  +  CH3  (k5  =  7.5  x  10~13  cm3  s-1)  (5) 

O  +  HCO  —  C02  +  H  (k6  =  1.0  x  10”10  cm3  s~‘)  (6) 

Note  that  the  energy  released  by  reaction  6,  AF/rxn  =  —467.8 
kJ  mol-1,  is  sufficient  to  populate  the  highly  excited  C02 
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Figure  6.  (a)  Log— log  plot  of  integrated  band  intensity  of  the  C02  V3 
emission  region  (1900—2000  cm-1, 0—50  us)  vs  the  laser  fluence  (upper 
curve)  and  SO2  density  (lower  curve),  (b)  Log-log  plot  of  the  chemical 
model  prediction  of  integrated  CO2  production  (0—50  ps)  vs  the  initial 
O-atom  density. 


vibrational  states  observed  in  the  data.  For  low  initial  O-atom 
concentrations  the  C02  chemical  production  rate  is  proportional 
to  [O]2  through  the  two-step  process  (5)— (6).  For  higher  O-atom 
concentrations,  the  net  HCO  production  rate  falls  to  zero  at  short 
times  following  the  photolysis  pulse  due  to  loss  via  reaction  6 
and  other  processes;  see  Figure  5  for  the  predicted  [HCO]  time 
dependence  for  [0]iniUai  =  6x  1014  cm-3.  In  the  latter  case  the 
C02  production  rate  is  linear,  assuming  the  HCO  concentration 
is  approximately  constant  over  the  time  integration  period  (0— 
50  ps).  The  chemical  model  reproduces  this  crossover  from 
quadratic  to  linear  behavior  in  [O], 

3.3.2.  Dependence  of  CO2  Emission  Intensity  on  [C2H4]. 
Figure  7a  shows  the  C02  emission  intensity,  integrated  between 
1900  and  2000  cm-1  and  0—25  ps,  plotted  as  a  function  of  the 
C2H4  concentration;  Figure  7b  shows  the  predicted  CO2  chemi- 
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Figure  7.  (a)  Log-log  plot  of  the  integrated  band  intensity  in  the 
CO2  V3  emission  region  (1900-2000  cm"1,  0-25  fx s)  vs  the  C2H4 
density,  (b)  Log-log  plot  of  the  chemical  model  prediction  of  integrated 
CO2  production  (0-25  n s)  vs  the  C2H4  density. 

cal  production  over  the  same  time  period.  The  C2H4  concentra¬ 
tion  range  in  Figure  7  is  about  1  —  10  times  that  of  [Ojjnidai-  In 
this  case  the  dependence  should  be  approximately  linear,  since 
the  CO2  production  rate  is  proportional  to  [C2H4]  at  early  times 
according  to  the  mechanism  given  by  reactions  5—6.  The 
average  slope  for  the  model  prediction  (0.94)  is  consistent  with 
this,  while  the  slope  for  the  data  is  marginally  lower,  at  0.73  ± 
0. 16.  The  smaller  value  for  the  slope  in  Figure  7a  could  be  due 
to  differential  quenching  of  the  CO2  emission  by  C2H4,  even 
during  the  0—25  [is  time  integration.  This  effect  is  discussed  in 
the  following  section. 

3.4.  Vibrational  Relaxation  of  CO2  by  C„H2 „.  For  variable 
amounts  of  C2H4,  the  integrated  radiance  in  the  1900-2000 
cm-1  region  reaches  a  maximum  at  delay  times  that  vary 
strongly  as  a  function  of  the  C2H4  concentration,  indicating  that 
C2H4  is  an  efficient  relaxer  of  the  vibrationally  excited  CO2. 
Figure  8  shows  the  signal  time  evolution  for  three  different  C2H4 
concentrations,  showing  how  the  maximum  depends  on  the  C2H4 
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Figure  8.  Semilog  plot  of  the  CO2  V3  integrated  emission  (1900— 
2000  cm"1)  vs  delay  time,  for  three  different  C2H4  concentrations.  The 
data  were  fit  to  eq  7  by  varying  the  quenching  rate  Q  and  the  scaling 
factor  C.  The  inset  shows  the  model  CO2  production  rate  Pit)  for  the 
three  C2H4  concentrations,  on  the  same  time  and  relative  intensity  scales 
as  the  data. 


concentration.  To  estimate  the  decay  rates,  each  curve  was  fit 
by  integrating  the  expression 


m 

dr 


=  CP(t)  -  QI(t) 


(7) 


In  eq  7,  /(f)  is  the  time-dependent  signal  intensity,  P(t)  is  the 
chemical  production  rate  from  the  chemical  kinetics  model  (in 
cm-3s_l),  Q  is  the  relaxation  rate  (in  s"1),  and  C  is  an  arbitrary 
constant  relating  the  radiant  intensity  to  the  number  density. 
Note  that  the  quenching  of  the  fluorescence  only  affects  the 
rising  portion  and  the  magnitude  of  the  population  curves.  The 
long-term  decay  rates  are  dictated  by  the  decay  of  the  feed  term 
P(t')  as  the  O  atoms  are  consumed  by  the  rate-determining 
reaction  5.  The  decay  of  the  feed  term  is  well-matched  by  the 
kinetic  model  with  no  adjustable  parameters,  as  shown  in  the 
inset  to  Figure  8. 

The  exponential  decay  constants  for  the  relaxation  of  highly 
vibrationally  excited  C02  as  a  function  of  ethene,  propene,  and 
n-butene  concentration  are  shown  in  Figure  9.  The  decay 
constants  were  obtained  by  fitting  the  1900—2000  cm-1 
integrated  emission  time  profile  to  eq  7  and  plotting  the  best- 
fit  disappearance  rates  Q  against  the  alkene  concentration.  The 
feed  functions  P(t)  for  propene  and  n-butene  were  assumed  to 
be  similar  to  that  for  ethene,  with  an  approximately  factor-of- 
five  shorter  time  scale  due  to  their  faster  reactions  with  O 
atoms.27  Propene  and  n-butene  were  found  to  give  rise  to  similar 
IR  spectral  profiles  in  reactions  with  O  atoms;  those  results  will 
be  detailed  in  a  future  publication.  As  discussed  above,  the 
1900—2000  cm-1  spectral  region  contains  emissions  from  highly 
excited  CO2  in  a  regime  where  cascade  from  higher  levels  is 
unlikely  to  present  a  major  effect.  For  higher  wavenumber  band- 
passes,  the  emission  occurs  from  lower  vibrational  levels  of  CO2 
that  are  more  prone  to  feed  from  higher  levels.  In  particular, 
nonexponential  decays  at  long  times  were  observed  when 
integrating  over  higher  wavenumber  ranges,  consistent  with  a 
more  complex  feed  mechanism  and  possibly  a  contribution  from 
radiating  CO  v. 
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Figure  9.  Best-fit  exponential  decay  constants  for  quenching  of  the 
integrated  C02  v3  emission  (1900-2000  cm'1)  as  a  function  of  the 
C2H4,  C3H6,  and  n-C4Hg  concentration. 


TABLE  3:  Rate  Constants  Icq  for  the  Quenching  of  C02 
Emission  (1900-2000  cm'1)  by  C2H4,  CdL,  and  n-C4Hg, 
Compared  to  Relative  Lennard- Jones  Collision  Frequencies 


kQ  (cm3  s'1) 

kQ(relative) 

fc0ii(L-J)(relative) 

C2H4 

(3.3  ±  0.8)  x  10'12 

1.00 

1.00 

c3h6 

(3.8  ±  0.5)  x  10'12 

1.15 

1.03 

/1-C4H8 

(4.2  ±  0.7)  x  10'12 

1.27 

1.12 

The  emitting  vibrational  levels  and  energies  can  be  ap¬ 
proximately  identified  by  using  the  known  anharmonic  and 
vibrational  mode  coupling  constants  of  CO2.  Emission  in  the 
1900—2000  cm-1  region  corresponds  to  roughly  15-19  quanta 
of  energy  in  the  v3  mode.  Of  course,  fewer  v3  quanta  are  excited 
if  the  exothermicity  is  partially  channeled  into  the  lower  energy 
V\  and  Vi  modes.  These  v3  levels  have  associated  term  energies 
of  33  000—41  000  cm-1,  as  much  as  93%  of  the  CO— O  bond 
dissociation  energy,  D0  =  43  980  cm-1. 

The  derived  rate  constants  for  the  relaxation  of  the  excited 
CO2  by  C2H4,  C3H6,  and  «-C4H8  are  given  in  Table  3.  The  values 
increase  only  modestly  with  increasing  carbon  chain  length,  with 
the  quenching  roughly  1%  efficient  based  on  a  gas  kinetic 
collision  rate.  The  increase  is  similar  to  the  estimated  12% 
increase  in  the  collision  rate  calculated  by  using  Lennard-Jones 
parameters  derived  from  experimental  second  virial  coef¬ 
ficients.35,36  The  alkenes  have  C=C  stretch  frequencies  (1623 — 
1647  cm-1)34  that  are  as  little  as  250  cm-1  lower  in  frequency 
than  the  high-energy  CO2  v3  vibrations,  suggesting  that  V— V 
transfer  could  play  a  role  in  the  relaxation  process.  Note  that 
C3He  and  n-C+Hg  have  IR-active  C=C  stretch  modes  with 
nonzero  transition  dipole  moments,  unlike  C2H4,  which  has  a 
center  of  symmetry.  However,  the  fact  that  the  relaxation  rate 
constants  are  not  significantly  larger  for  C3H6  and  n-GrHg 
implies  that  V— V  transfer  induced  by  long-range  forces  does 
not  make  a  significant  contribution  to  the  quenching  mechanism. 
Relaxation  by  SO2  was  found  to  be  very  inefficient  (<10~15 


cm  3  s  '),  consistent  with  the  lack  of  a  vibrational  mode 
frequency  near  those  of  highly  excited  CO2. 

4.  Summary 

A  new  experiment  has  been  initiated  to  measure  and 
characterize  the  IR  emissions  resulting  from  the  interactions  of 
O  atoms  with  hydrocarbons.  In  an  initial  study,  O  atoms  were 
prepared  by  photolyzing  S02,  then  reacted  with  C2H4.  The 
nascent  fast  O  atoms  are  largely  thermalized  through  collisions 
with  Ar  bath  gas  prior  to  reaction,  resulting  in  O  +  C2H4  primary 
reaction  and  related  secondary  chemistry  proceeding  under 
thermal  conditions.  The  resultant  IR  chemiluminescence  signal 
was  detected  by  using  a  step-scan  FTIR  spectrometer  for 
simultaneous  time-  and  frequency-dependent  data  collection. 
Emission  has  been  observed  from  CO,  CO2,  HCO,  H2CO,  and 
other  product  species.  Integrated  radiance  in  the  CO2  v3  emission 
region  has  been  correlated  with  reactant  species  partial  pressures 
and  laser  fluence.  A  chemical  kinetics  model  has  been  used  to 
identify  secondary  chemistry  responsible  for  the  appearance  of 
chemiluminescent  CO  and  C02  and  to  characterize  the  partial 
pressure  and  fluence  dependences.  Rate  constants  for  the 
relaxation  of  highly  vibrationally  excited  CO2  by  C2H4,  C3H6, 
and  n-C4Hg  have  been  determined,  the  values  increasing  only 
modestly  with  increasing  alkene  chain  length. 
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